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Abstract The climatic preferences of Pelias barani, a rare Pontic endemic viper, are analysed and a new locality record 
is reported. According to species distribution modeling, with the average test AUC was was 0.904 + 0.068, bioclimatic 
variables such as Precipitation of driest month (45.5%), Temperature seasonality (18.9%), Precipitation seasonality 
(17.9%) and Maximum temperature of warmest period (14.7%) appear to have the most useful information on 
geographic distribution of Pelias barani. Distribution models of Pelias barani under current climatic conditions showed 
better adaptation to the northwest and northeast part of Turkey. 
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1. Introduction 


The family Viperidae includes three subfamilies, thirty- 
eight genera (Pyron et al., 2013) and 328 species (from 
The Reptile Database 2014). The genus Vipera (Laurenti, 
1768) includes 32 species inhabiting Asia Minor, Eurasia, 
Northern Africa and the Middle East (Ananjeva et al., 
2006; from The Reptile Database 2014). Currently, 14- 
15 species belonging to the family Viperidae have been 
recognized in Turkey (Tuniyev et al, 2012, Göçmen et al., 
2014). Taxonomically, the genus Pelias was recognized 
first as a subgenus of genus Vipera (Venchi and Sindaco, 
2006), and many authors used “Pelias” as subgenus 
(Ananjeva et al., 2006; Venchi and Sindaco, 2006; Afsar 
and Afsar, 2009; Avcı et al., 2010), but later according 
to a catalogue of living and extinct species in recent 
years (Wallach et al., 2014), Pelias was recognized as 
a genus in the family Viperidae that includes P. altaica, 
P. anatolica, P. barani, P. berus, P. darevskii, P. dinniki, 
P. ebneri, P. eriwanensis, P. kaznakovi, P. lotievi, P. 
magnifica, P. nikolskii, P. olguni, P. orlovi, P. pontica, 
P. renardi, P. sachalinensis, P. seoanei, and P. ursinii 
species (Wallach et al., 2014). Accordingly, I use Pelias 
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as genus in this study. Pelias barani (Baran’s Adder) 
is one of the Anatolian vipers with little information 
available about its biology (Kumlutas et al., 2012). P. 
barani was first described from its type locality in 60 
km N of Adapazarı at 400 m a.s.l. (40°50' N, 30°25' E) 
in the northwestern Anatolia with a female specimen 
(Böhme and Joger, 1983). Later, new locality records 
were reported from the northwestern and the northeastern 
Anatolia (Baran et al., 1997; Franzen and Heckes, 2000; 
Baran et al., 2001; Avcı et al., 2004; Baran et al., 2005; 
Kumlutas et al., 2012) (Table 1). The taxonomic status of 
P. barani was evaluated with several studies in respect to 
phylogenetic data. Firstly, Joger et al. (1997) studied the 
phylogenetic position of P. barani and P. nikolskii within 
the P. berus complex using morphological and hemipenial 
comparisons, and a partial sequence of the mitochondrial 
cyt b gene. According to mitochondrial DNA information, 
they found that P. barani might be a subspecies of P. 
berus, but they showed that it was a different species 
using both morphological and hemipenial data. Later, 
Joger et al. (2003) identified five groups including P. 
barani, and genus Pelias that is consisted basically of 
two lineages covering P. berus and P. ursinii-kaznakowi 
complex. Secondly, Kalyabina-Hauf et al. (2004) studied 
systematic phylogeny of P. berus complex, and they found 
that P. bosniensis and P. barani were morphologically 
different from P. bersus and are suspected to be a 
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Table 1 Morphometric measurements and counts of known specimens based on literature information together with new locality. SSO: Scales between supraoculars, EUL: Scale rows between 
eyes and upper labials, DS: longitudinal rows of dorsal scales. References: 1) Böhme and Joger (1983), 2) Baran et al. (1997), 3) Franzen and Heckes (2000), 4) Baran et al. (2001), 5) Avci et al. 


(2004), 6) Baran et al. (2005), 7) Kumlutas et al. (2012), 8) This study 
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different species of them. In addition, they reported that 
if P. barani would be recognized a separate species, P. 
barani was supported as paraphyletic, with P. berus. 
Finally, Garrigues et al. (2005) indicated that molecular 
phylogeny of genus Vipera, which has three main 
branches of the European clade, were distinguished the 
Pelias group including P. barani, Vipera ammodytes, and 
Vipera aspis. 

Species distribution modeling creates an expected 
distribution map of a species based on climatic and others 
environmental conditions known presence localities 
(Kozak et al., 2008). It is also useful for conservation 
strategies and selection of protected areas (Gül, 2013) 
as well as the relationship between divergence and 
speciation mechanisms, potential geographic distributions 
of species on speciation (Graham et al., 2004). 

In the text the main goal is to define the climatic 
patterns and factors affecting species distribution. 
Therefore, I ran species distribution modeling by 
Maxent (Philips et al., 2006) because knowledge about 
the geographical distribution of P. barani is crucial for 
conservation and spatial planning. 


2. Materials and Methods 


2.1 Studied Species A female Pelias barani specimen 
was found in village Büyükköy of Çayeli, Rize, at 529 
m above sea level, on October 18, 2014 by Rahime 
USTABAS from Twelfth Region Headship in Rize under 
Ministry of Forest and Water Affairs (Figure 1). After 
analysis of morphological features of the specimen (Table 
1) at the Zoology Research Laboratory, Recep Tayyip 
Erdogan University, it was safely released back into its 
natural habitat. 


2.2 Environmental Variables 19 bioclimatic variables 
were downloaded from global climate layers (available at 
www.worldclim.org) in the highest resolution at 30 arc- 
seconds (~1 km) under current conditions (~1950—2000) 
(Hijmans et al., 2005). The global climate layers were 
clipped to the borders of Turkey using Extract by Mask in 
ArcGis v. 10.1. Many of the bioclimatic variables are very 
similar to each other. Therefore, all bioclimatic variables 
were examined for Pearson correlation coefficient (0.75 
< r <-0.75) using ENMTools 1.3 (Warren et al., 2010) 
and redundant variables were excluded. Six bioclimatic 
variables were eventually selected for the model (Table 2). 

Species distribution modeling was performed using 
Maxent software v. 3.3.3e (Philips et al., 2006). In order 
to develop species distribution modeling, 10 presence 
localities were used based on the new locality record 
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Figure 1 Overview of female Pelias barani from Büyükköy, 


and literature data (Figure 2, Table 3). Seven test points 


similar to presence localities were created by using 
Geospatial Modeling Environment version 0.7.2.0 (Beyer, 


2012). The final model was composed of the average the 


3.1 Morphology of the new specimen of Pelias barani 
The number of ventral scales of the specimen caught 
from Çayeli is 142 and the number of subcaudal scales is 
29-30. The scales on longitudinal rows of dorsal at mid- 
body were 21. The specimen had two apicals in contact 
with rostral and also had two canthals on each side of 
the head. Loreal scales between the preocular and the 
postnasal were 5-5. Circumocular scales around the eye 
were 11—12 in the specimen. Upper and lower labials of 
the specimen were 9—9 and 11—12, respectively. Gular 
scales that are in contact with the first ventral scales were 
4—4. The specimen had a total length of 670 mm (head 
and body length 600 mm; tail length 70 mm) (Table 1). 

3.2 Color patterns of Pelias barani The head of the 
female specimen captured in Cayeli is large (Figure 
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Kilometers 


Figure 2 This picture is a representation of the model for Pelias barani based on both literature data and new locality record. Warmer colors 


show areas with better predicted conditions. 


Table 3 Localities information of Pelias barani based on literatures data. 


Species Latitude Longitude Localities References 

P. barani 41.172167 30.227777 60 km N of Adapazari Bohme and Joger, 1983 
P. barani 41.0462273 37.0433088 İkizce, Ordu Avcı et al., 2004 

P. barani 40.584711 40.410233 Arpagözü Yaylasi, Çaykara, Trabzon Baran et al., 2001 

P. barani 40.6090956 30.5564713 Göktepe, Geyve, Adapazarı Baran et al., 2005 

P. barani 41.4129954 31.7502575 Kozlu, Zonguldak Kumlutas et al., 2013 

P. barani 40.7397124 38.4615057 Dereli, Giresun Franzen and Heckes, 2000 
P. barani 41.105277 41.044452 Duygulu köyü, Çamlıhemşin, Ardeşen Baran et al., 1997 

P. barani 41.0487107 41.007611 Fırtına valley, Ardeşen, Rize Franzen and Heckes, 2000 
P. barani 41.4699826 41.8979419 The northwest of Artvin Baran et al., 2001 

P. barani 40.994275 40.691457 Büyükköy, Çayeli, Rize This study 


3A). The dorsal color pattern is almost black in hue. 
The ventral color includes many different shades of 
brown, sometimes darkish, or whitish brown and also 
ground color of the ventral side is lighter than dorsal side 
(Figure 3B). This color variation continues across both 
upper labials and lower labials on each side of head. The 
specimen has a different color of ventral surface of the 
end of the tail that is orange (Figure 3C). In addition, 
references, which are listed in Table 1, indicated that the 
body color of all known specimens were totally black; 
however, the ground color of the second specimen of 
reference 4 was greyish brown with a blackish zig-zag 
band. 


3.3 Species distribution modeling Based on current 
climatic conditions distribution model of Pelias barani 
showed that some parts of northwest and northeast of 
Turkey was the areas represented with better predicted 
conditions for Pelias barani’s habitat (Figure 2). 
As a result of estimates of relative contributions of 
the environmental variables with Maxent, the most 
environmental variables with highest gain that explains 
more than 10% of the presence of P. barani were 
Bio-14 (Precipitation of driest month, 45.5%), Bio-4 


(Temperature seasonality, 18.9%), Bio-15 (Precipitation 
seasonality, 17.9%) and Bio-5 (Maximum temperature of 
warmest period, 14.7%). Other variables had a percent 
contribution less than 5%. The average test AUC (the 
area under the receiver operating characteristic curve) 
value of the distribution model of P. barani was 0.904 


+ 0.068 and 10 percentile training presence logistic 
threshold was 0.3632. Predicted suitability was highest in 
northwest and northeast regions of Turkey with extreme 
precipitation conditions that influence potential range of 
P. barani during the year, but with a lower temperature 
variation throughout the year, low variation in total 
monthly precipitation throughout the year and low warm 
temperature anomalies throughout the year (Figure 4). 


4. Discussion 


There are two-distribution patterns for P. barani in 
Turkey, including the Black Sea Region (Zonguldak, 
Ordu, Giresun, Trabzon, Rize and Artvin localities) and 
the Marmara Region (Adapazarı localities). The Black Sea 
Region has the greatest amount of precipitation in Turkey. 
Its summer is warm and humid and the winter is cool and 
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Figure 3 Dorsal (A), and ventral (B) side of typical head pattern, 
and color of the ventral surface at the end of tail (C) in female 
specimen from Büyükköy, Çayeli. 


damp (Borysova et al., 2005). Marmara Region consists 
generally of three zones (istanbul and Yildiz Mountain to 
the north, izmir and Ankara to the south, and Sakarya). 
Therefore, it is the region that has the greatest diversity 
of climate variable. It has Mediterranean climate that is 
warm to hot, dry summers and mild to cool, wet winters 
on the Agean Sea coast and the south Marmara Sea coast, 
and has an oceanic climate that is warm summers and 
cool winters on the Black Sea coast (Sertel and Ormeci, 
2011). Adapazari localities found in the Sakarya zone 
of P. barani are similar to the Black Sea Region with 
regard to climatic features (Sirdas, 2005). The geographic 
distribution of P. barani that is limited with the north 
coastal areas of Turkey indicates that this species is not 
resistant to a wide variety of conditions (eurioic) (Figure 
2). P. barani has actually been detected in particular 
regions under very rainy environmental conditions on the 
Anatolian range. I found that the model, which explains 
the current distribution pattern of P. barani, is based on 
the precipitation of driest period (Bio-14), Temperature 
seasonality (Bio-4), Precipitation seasonality (Bio-15) 
and Maximum temperature of warmest period (Bio-5) of 
northwest and northeast of Turkey. Summer temperature 
and precipitation are probably a limitation factor of P. 
barani’s dispersion (Figure 2). Therefore, P. barani 
inhabits towards to the Black Sea region of Turkey. 
In fact, the geographical structure of Anatolia plays a 
significant role for a rich biodiversity because the Anatolia 


generates from multiple glacial refugia (Bilgin, 2011; 
Tarkhnishvili et al., 2012). This suggests that occurrence 
of the genetic diversity correlates with differentiation in 
the ecology of local refugia within Anatolia (Giil, 2013). 
The important mountain chains, such as Toros Mountains, 
North Anatolian Mountains, the Anatolian Diagonal and 
Western Anatolian Mountains, are a barrier affecting the 
distribution of many amphibians (Gül, 2013; Ozdemir 
et al., 2014) and reptiles (Kapli et al., 2013; Sindaco 
et al., 2014). Similarly, Brito et al. (2008) showed that 
distribution of the viperid snakes (Vipera latastei and V. 
monticola) were influenced by the precipitation of 
driest period for Eastern Iberia, Algeria, and Rif and 
Middle Atlas, and High Atlas except Western Iberia. 
In biogeographical scenario of Vipera seoanei, another 
study represents that the origin of genus Pelias can be 
expanded from the north Black Sea region toward Europe 
between the late Miocene and early Pliocene, and range 
expansions of the genus Pelias were affected by cooler 
temperatures of Pleistocene while warm temperatures of 
middle Pliocene prevailed its expansion (Martinez-Freiria 
et al. 2015). In fact, species of the genus Pelias can be 
better adapted to colder climates (Garrigues et al., 2005). 
In addition, several studies support that climatic factors 
are responsible for geographical distribution of many 
reptile species (Santos et al., 2006; Brito et al., 2008; 
Brito et al., 2011; Sow et al., 2014; Martinez-Freiria et 
al., 2015). 

Consequently, the P. barani specimen with new 
locality in this study has similarities to specimens in 
the literature in terms of morphometric measurements 
(Table 1). Since P. barani is a rare species, new locality 
records are very important in order to understand the 
distributional ranges of the species. In addition, this study 
represents that climatic conditions restricts its potential 
distribution in Turkey, and also it has been suggested 
that local precipitation patterns was important for the 
distribution model of the species. 

According to IUCN Red List, P. barani (Vipera barani 
in IUCN) is listed as near threatened (NT) and it is 
known as pontic endemic (Venchi and Sindaco, 2006). 
Population trend of P. barani has been decreased because 
of habitat loss, and persecution. Similarly, I observed 
habitat loss caused by destruction of forest through com- 
mercial harvest and construction of roads in Büyükköy 
locality. In addition, it is killed like many snakes when 
P. barani is seen by the local people. I also realized that 
many amphibians and reptiles were being killed on roads, 
and thought that P. barani might be killed like others. 
Therefore, habitat destruction, road deaths and deaths 
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Figure 4 Response curves of the four bioclimatic variables that limit the distribution of P. barani. (A) Precipitation of driest period (Bio14), (B) 
Temperature seasonality (Bio4), (C) Precipitation seasonality (Bio15) and (D) Maximum temperature of warmest period (BioS). 


caused by local human activities is a major threat to 
survival of P. barani and significant cause of extinction 
of P. barani throughout Black Sea coast. As a result of all 
these threats, it might be candidate to enter Vulnerable 
(VU) category under the red list category in the near 
future. Literature information about its distribution and 
ecology are very scarce. Therefore, it’s a new locality 
record and climatic preferences might be useful for 
conservation of this pontic endemic viper. 
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